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ABSTRACT 

This paper describes the algorithms of the Photometer Performance Assessment (PPA) software component in the 
science data processing pipeline of the Kepler mission. The PPA performs two tasks: One is to analyze the health and 
performance of the Kepler photometer based on the long cadence science data down-linked via Ka band approximately 
every 30 days. The second is to determine the attitude of the Kepler spacecraft with high precision at each long cadence. 
The PPA component is demonstrated to work effectively with the Kepler flight data. 
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1. INTRODUCTION 

The Kepler mission is the first NASA mission capable of detecting Earth-size planets in the habitable zone of solar-like 
stars 1 ' 2 . The mission consists of a spacecraft with a single scientific instrument, a photometer, which is used to 
continuously and simultaneously measure the brightness variation of target stars to look for sequences of planetary 
transits. The Kepler spacecraft was launched into an Earth-trailing heliocentric orbit from Cape Canaveral, Florida on 
March 6, 2009. The baseline mission is planned for 3.5 years of operation following a 65-day commissioning period. 

1.1 Photometer of the Kepler spacecraft 

The photometer consists of a classic Schmidt telescope design with a 95 cm aperture and a field-of-view (FOV) of 
-115 square degrees. The focal plane array consists of 42 Charge-Coupled Devices (CCDs). Each CCD has 1024x2200 
pixels and two outputs with the serial channel on the long edge. Two CCDs are mounted together with a field flattener 
lens to form a nearly square module. As a result, there are a total of 21 modules, 84 outputs, and 96 megapixels on the 
Kepler photometer focal plane. 

Each CCD integrates for 6 seconds and is read out in 0.5 seconds. The data from each pixel are co-added onboard the 
Kepler spacecraft for 30 minutes. Each of these 30-minute intervals is known as a long cadence, and the 30-minute co- 
added data as long cadence data. Using an onboard target table, data from the pixels of interest for each of -156,000 
target stars, along with collateral data used for calibration are extracted, compressed, and stored onboard. 

During scientific observation of the entire mission, the attitude of the Kepler spacecraft is three-axis stabilized to keep 
the photometer FOV at a fixed region of the sky. The downsampled data of a sparse set of selected targets, which are 
used to monitor the photometer’s health and status, are downloaded to the ground semi-weekly with an omni-directional 
X-band (8 GFlz) antenna, which avoids interruption of scientific observation. Approximately every 30 days, scientific 
observation is halted and the spacecraft is maneuvered to point a Ka-band (32 GFlz) high-gain antenna toward the Earth 
to download all the data in storage to the ground. In order to maintain the solar array toward the Sun, the spacecraft is 
rotated 90 degrees about the optical axis of the photometer every three months. The CCD layout on the focal plane is 
designed to be four-fold symmetric and CCD orientations are set up so that the readout row and column directions 
remain the same with each rotation (except for the center module). 
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1.2 Photometer performance assessment software component 

Photometer Performance Assessment (PPA) is one of the software components of the science data processing pipeline of 
the Kepler mission' 4 . One task of PPA is to analyze the health and performance of the Kepler photometer. The 
instalment metrics, either determined within the PPA component or provided as the outputs of the Calibration (CAL) 5 or 
Photometric Analysis (PA) 6 components, are analyzed with a track and trend algorithm where the mean and uncertainty 
are determined, out-of-bound values are identified and the trend is predicted. Another task of PPA is to determine the 
attitude of the Kepler spacecraft from the measured centroid position of the images of bright target stars. To compress 
the data to be processed, two-dimensional motion polynomials are built in the PA component to map the right ascension 
and declination of a target star to the row and column of the centroid of the star image for each output at each long 
cadence. The motion polynomials for all outputs are provided as inputs to PPA, which uses them to determine the 
attitude with high precision. 

Compared to the Photometer Data Quality (PDQ) component 7 which generates a report of the photometer’s health and 
status from a subset of data downloaded via X-band semi-weekly, the PPA component analyzes the photometer’s health 
and performance and determines the spacecraft’s attitude from all the long cadence data downloaded via Ka-band every 
30 days. 

This paper is organized as follows: First, the Kepler instrument metrics are discussed, and the track and trend algorithm 
and the photometer performance report of the Kepler flight data are presented; Second, the data compression and 
smoothing with motion polynomials are described, the attitude determination algorithm with nonlinear least squares fit is 
described, and the attitude solutions of the Kepler flight data are demonstrated; And last, conclusions are presented. 


2. PHOTOMETER PERFORMANCE ANALYSIS 

2.1 Kepler instrument metrics 

In the Kepler science data processing pipeline, data are processed per output 8, 9 . In the PPA component the health and 
performance of the Kepler photometer are assessed by analyzing metric time series such as general instalment metrics, 
cosmic ray metrics, Combined Differential Photometric Precision (CDPP) metrics, etc. for each of 84 outputs and two 
additional metric time series - theoretical and achieved compression efficiencies for the fill focal plane. 

Table 1 lists general instalment metrics for each output and provides the unit and a brief description for each metric. The 
software component, in which each metric is determined, is shown in parenthesis in the description. For example, the 
metrics of black level, smear level and dark current are determined in CAL 5 and the metrics of brightness and encircled 
energy are determined in PA 6 . 

The background level is one of the general instrument metrics determined in PPA. The background values of an 11x11 
grid of pseudo targets are determined for each output by evaluating the background polynomial provided by PA. The 
background level metric of the output is calculated as the weighted mean of the background values of the pseudo targets, 
where the weight is inversely proportional to the square of the uncertainty of each value. Figure 1 shows the background 
level metrics across the focal plane of the flight data on the Modified Julian Day (MJD) 55023, or July 11, 2009. The 
background level increases gradually from the lower left to the upper right comer of the figure, since the upper right 
corner is closer to the galaxy equator (where more dim stars are located), as a result of the spacecraft’s quarterly roll. 

In PPA, theoretical and achieved compression efficiencies of the full focal plane are determined and analyzed as metrics 
of efficiency of the onboard data compression algorithm, as well as metrics of the margin of the memory usage of the 
onboard computer. Compression efficiencies of the full focal plane are calculated from the corresponding metrics of 
each output, provided by the CAL component, with the following equation: 

84 

2>< -CE, 

CF - -is! 

2 >< 

i = 1 


(i) 



where s t is the number of code symbols of output i (7=1,. . .84), CE, is the theoretical or achieved compression efficiency 
of the output, and CE FFP is the corresponding compression efficiency of the full focal plane. 


Table 1. General instalment metrics for each of 84 outputs. 


Metrics 

Unit 

Description 

black level 

ADU/read 

Median of black level (CAL) 

smear level 

electron/cadence 

Median of smear level (CAL) 

dark current 

electron/second 

Median of dark current (CAL) 

brightness 

dimensionless 

Mean of the ratio of received flux to the expected flux of the target star 
(PA) 

encircled energy 

pixel 

Mean of the radius of the circle to hold 95% of the energy of the target 
star (PA) 

background level 

electron/cadence 

Median of background level (PPA) 

centroid mean row 

pixel 

Mean of difference in the row between the centroid of the image of the 
target star and its predicted position from nominal attitude (PPA) 

centroid mean column 

pixel 

Mean of difference in the column between the centroid of the image of 
the target star and its predicted position from nominal attitude (PPA) 

plate scale 

arcsecond/pixel 

Mean of the ratio of angular distance between stars and the linear 
distance between their images on the focal plane (PPA) 

theoretical 

compression efficiency 

bit/pixel/cadence 

Theoretical compression efficiency of the onboard data compression 
algorithm (CAL) 

achieved 

compression efficiency 

bit/pixel/cadence 

Actual achieved compression efficiency of the onboard data 
compression algorithm (CAL) 


x 10 6 



Figure 1. Background level across the focal plane of the flight data for MJD 55023 (unit: electron/cadence). 



Figure 2 shows achieved compression efficiency of the full focal plane of flight data for second quarter of Kepler science 
data collection, ranging from MJD 55002 (June 20, 2009) to MJD 55091 (September 17, 2009). Spikes are observed in 
the achieved compression efficiency after events of attitude variations such as a safe mode, an attitude maneuver, or an 
attitude tweak that lead to a change in the star images. As a result, redundancy of raw data decreases and size of 
compressed data increases. 
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Figure 2. Achieved compression efficiency of the full focal plane for second quarter flight data 

Cosmic ray metrics are a set of time series for the hit rate of cosmic ray and the mean and variance of the distribution of 
cosmic ray energy for a specified group of pixels on the CCD. For each output, CAL generates cosmic ray metrics of 
black pixels, masked smear pixels, and virtual smear pixels 3 , and PA generates cosmic ray metrics of target star pixels 
and background pixels 6 . 

CDPP is a metric for the effective white noise standard deviation in a time interval for the flux time series of the target 
stars. It determines the signal-to-noise ratio of the transit with a given depth and duration of that interval. In PPA, target 
stars are grouped in seven bins according to their magnitude, from Magnitude 9 to Magnitude 15, for each output. The 
expected and measured CDPP metrics are determined for each magnitude bin from the flux value and uncertainty time 
series, generated by the Pre-search Data Conditioning (PDC) component 10 , and the CDPP time series, generated by the 
Transiting Planet Search (TPS) component 11 , respectively. 

2.2 Track and trend algorithm 

The time series of the instrument metrics are analyzed with the track and trend algorithm, which identifies metric 
samples beyond the fixed/adaptive upper and lower bounds, and predicts whether the metric crosses the fixed/adaptive 
bounds within a horizon time beyond the end of the metric time series. Note that the fixed upper and lower bounds of the 
metrics and the horizon time for prediction are pre-defined input parameters of PPA. 

Adaptive upper and lower bounds are determined from estimates of time-varying mean and noise/variation levels of 
metric time series. For a given metric time series {dt, k= 1 , 2, . . . }, the mean value, variance and uncertainty are estimated 
with a first order lowpass filtering algorithm as follows: 
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where <4 is the metric sample at the k - th long cadence, m l: , v/ f and lit are the estimated mean value, variance and 
uncertainty at the same long cadence respectively. /?,„ and /?,,, ranging between 0 and 1, are filtering coefficients for the 
mean value and variance respectively. A filtering coefficient closer to 0 suppresses the noise more effectively at the cost 
of sluggishness to follow the time variation of the metric time series. Based on the estimates, the adaptive upper and 
lower bounds are defined by the time series {nik+Uk) and {ffirW/t} respectively. 

The trend of metric time series is determined by selecting metric samples in a pre-defined interval (for example, five 
days) at the end of the time series and then fitting the selected data with a linear model. It is straight-forward to 
determine when the trend line specified by the linear model crosses the fixed and adaptive bounds. 

An alert is generated, attached with an alert level, when a metric sample is identified to be beyond fixed or adaptive 
bounds or when there is a bound-crossing event within the horizon time beyond the end of a metric time series. For 
example, samples beyond the fixed and adaptive bounds are marked red and yellow respectively, and the alert level is 
marked green when no samples are identified beyond the fixed or adaptive bounds, 

Figure 3 shows a diagram of the track and trend algorithm and Figure 4 illustrates the track and trend of achieved 
compression efficiency of the full focal plane for second quarter flight data. In Figure 4, the metric samples are plotted 
with a solid line and dot markers, the adaptive upper and lower bounds are plotted with dash lines, and the trend is 
plotted with a dash-dot line. The estimated uncertainty increases to indicate higher noise/variation levels when the spiky 
metric samples are encountered, and the band defined by the adaptive upper and lower bounds becomes wider 
accordingly. When the estimated uncertainty converges to the normal noise/variation level, the band shrinks, showing 
adaptation to the time-varying feature of the mean and noise variation level of the metric time series. 



Figure 3. The track and trend algorithm of instrument metrics. 
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Figure 4. Track and trend of achieved compression efficiency of the full focal plane for second quarter flight data 

2.3 Photometer performance reports of flight data 

Photometer performance reports are generated by PPA for each output and for the full focal plane. Figures 5a and 5b 
show the general instrument metrics dashboard of second quarter flight data for module 6/output 3 and the full focal 
plane respectively. As seen in Figure 5a, each of the general instrument metrics corresponds to a rectangle color-coded 
with the track alert level. Figure 5b is an aggregation of the dashboard of all 84 outputs on the focal plane, where module 
6/output 3, is the first output of the third row. Each small rectangle in the square of an output corresponds to a metric, as 
labeled in Figure 5a. The two squares in the upper right corner of Figure 5b correspond to the metrics of theoretical and 
achieved compression efficiencies of the full focal plane. It is seen that the track alert levels of most metrics for the 
outputs across the focal plane are green, demonstrating that the Kepler photometer is in good health and performs well in 
the second quarter. The red alert level of achieved compression efficiency of some outputs is due to the spiky samples 
following the attitude variations, as shown in Figure 2. The in-flight performance of the Kepler instrument is discussed 
in detail by Caldwell, D.A., et al. 12, 13 


3. ATTITUDE DETERMINATION WITH HIGH PRECISION 

3.1 Measurement data compression and smoothing 

In PPA, the attitude of the Kepler spacecraft is determined from the measured centroid position of the images of bright 
target stars. Approximately 200 bright target stars are selected for each output. Based on the measurements of -16,000 
star images across the focal plane, the attitude solution of the Kepler spacecraft is determined with high precision at each 
30-minute long cadence. However, the size of the measurement data for a month or a quarter, which is the length of data 
processed by the Kepler science data processing pipeline, is so large that it is difficult to process the data directly to 
determine the attitude. As a solution to the problem, two motion polynomials, one for row and one for column, are 
generated in the PA component for each output and each long cadence. By fitting the measured rows or columns of the 
centroids of the images of ~200 bright target stars from an output into a two-dimensional polynomial of the 
corresponding right ascensions and declinations of the stars, the motion polynomial builds an empirical mapping from 
the position of a target in the sky to the measured position of the target on the focal plane of the Kepler photometer. 


(a) 



Figure 5. General instrument metrics dashboard for second quarter flight data, (a) module 6/output 3; (b) full 

focal plane. 


In PPA, a set of pseudo stars is selected for each output first, and the positions of the selected pseudo stars in the sky are 
determined; then, the “measured” rows and columns of the pseudo stars, together with the covariance matrices, which 
are calculated from the row and column motion polynomials, are used as input measurements for the attitude 
determination. Considering the order of the motion polynomials is limited (usually <4), the size of the “measurement” 
data of the pseudo stars, which is input to the attitude determination function in PPA, is reduced significantly and the 
noise in the measurements of original target stars is suppressed effectively as well. 

For each output, pseudo stars are selected from a 3 X 3 grid. The requirement that input measurements should be 
uncorrelated imposes constraints on the selection of pseudo stars. The correlation of rows or columns of a set of pseudo 
stars, determined by evaluating the two-dimensional motion polynomial, is related to the order of the motion polynomial 
and the geometry of the pseudo stars. Figure 6 illustrates the selection of pseudo stars of an output from a 3 X 3 grid in 
case of (a) N mp = 0, (b) N mp = 1, (c) N mp = 2, and (d) N mp >= 3 respectively, where N mp is the order of the motion polynomial. 

The selected pseudo stars for each of 84 outputs are fixed in the sky. Their positions, in right ascension and declination, 
are determined from the nominal attitude of the Kepler spacecraft, which will be discussed later. 

3.2 Attitude determination with nonlinear least squares fit 

The attitude of the Kepler spacecraft is described by the right ascension and declination of the optical axis of the 
photometer, denoted as a s and 8 S respectively, and a rotation angle along the optical axis, denoted as (j) s . A 3-2-1 

transformation with Euler angles of OC s , - 8 S and (f> s rotates the inertial equatorial coordinate system, with the X axis 
pointing to the vernal equinox, to the spacecraft body fixed photometer focal plane coordinate system; thus, pointing the 
photometer FOV to a fixed region of the sky. The nominal values of CC S and 8 S are 290.66667° and 44.5° for the 

Kepler science mode, which are constant for the entire mission. The rotation angle (f) s varies by multiples of 90° 

corresponding to quarterly spacecraft roll maneuvers, with a nominal value of 123° in the quarter from June 20, 2009 
(MJD 55002) to September 17, 2009 (MJD 55091). 









(C) (d) 

Figure 6. Selection of pseudo stars of an output from a 3x3 grid in case of (a) N mp =0, (b) N mp = 1 , (c) N mp = 2, and 
(d) N mp >=3 respectively, where N mp is the order of the motion polynomial. 


Let f/ l j and d t • denote the right ascension and declination of the /- 1 h selected pseudo star of output i, where i=\, ... 
84, y'=l, ..., J h Jj is the total number of selected pseudo stars of output i. As discussed in Section 3.1, the “measured” row 
and column of the centroid of the image of the pseudo star at the k-\h long cadence, denoted as r i . k and 

C . j k respectively, are determined by evaluating the row and column motion polynomials at <2 ; . and 8 i . , as expressed 
in the equations below: 


r i,j,k= MP row ,i,k( a i,j’ 5 ij) 


(3a) 


c ij,k= MP coi,i,k( a iJ’ S iJ 
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where MP l k (-,-) and MP col ^ k { •) are two-dimensional motion polynomials of the row and column for output i at 
the k- th long cadence. The covariance matrices of the measurements r j . k and C ( - , k (j= 1, . . .J,), denoted as W row t k and 
W co , t k , are determined from the covariance matrices of the coefficients of MP mw . k (•, •) and MP col . k (•, •) . 

The “predicted” row and column of the centroid of the image of the pseudo star at the k- th long cadence, denoted as 
r { j k and C ( ■ k respectively, are functions of the right ascension <2, . and declination . of the pseudo star, and the 

attitude parameters at the /c-th long cadence, denoted as a s k , S s k and (j) s k , as following: 


K j,k = MD r OW ,i ( a u , , a Sr k , S sk (/> sk ) 

P,j,k = MDcotJ ( a i,j ’ ^i,j ’ a s,k ’ ^s.k^S.k ) 
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where MD row . and MD col . denote the models of the spacecraft attitude and the photometer optics, including models 
of stellar velocity aberration, photometer optical distortion and the focal plane geometry etc. 14 ' L \ 




The attitude solution at the k - th long cadence is determined as the attitude parameters OC s k , S s k and (f) s k which 
minimize the nonlinear least squares function as follows: 
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Figure 7 shows a diagram of the attitude determination algorithm. As depicted in the diagram, the attitude determination 
algorithm is implemented in the following steps. First, a set of pseudo target stars is selected for each output based on the 
orders of the corresponding motion polynomials. The right ascension and declination of the pseudo star are determined 
from the nominal attitude of the Kepler spacecraft. Second, the “measured” rows and columns of the centroids of the 
pseudo stars of each of the 84 outputs at a 30-minute long cadence are calculated by evaluating the corresponding row 
and column motion polynomials, as shown in Equation (3). The covariance matrices of the measurements of the row and 
column are determined as well from the covariance matrices of the coefficients of the motion polynomials. Then, the 
outliers in the “measured” rows and columns are identified and removed with robust regression, and the attitude solution 
is determined at the long cadence as the nonlinear least-squares fit of the cleaned “measured” data of the pseudo stars 
across the focal plane, as shown in Equation (5). Finally, a metric of the maximum position offset of star images on the 
focal plane between the attitude solution and the nominal attitude is determined to see how well images of target stars 
are located within corresponding apertures defined by the target table. 

3.3 Attitude solutions of flight data 

Figure 8 shows the difference in declination between the nominal attitude and the attitude solution for second quarter 
flight data. The timestamps of the figure are elapsed days referenced to MJD 55002.0, or OhOmOs June 20, 2009. Several 
events of the Kepler spacecraft are observed in this figure. The spacecraft goes into safe mode on day 12 and comes back 
to normal science orientation on day 14. A transient in the attitude is observed following the safe mode, mainly due to 
temperature variations. On days 54 and 86, the spacecraft goes into coarse pointing mode. During days 61 and 62, the 
spacecraft is maneuvered to point its high-gain Ka-band antenna toward the earth and science data is downloaded to the 
ground. A pointing drift is observed in the measured declination relative to the nominal science attitude 16 . In order to 
keep all selected target stars within the corresponding apertures defined by the target table, the drift is corrected by 
attitude tweaks on days 0, 31, and 77. As expected, jumps in the spacecraft attitude parameters are observed during the 
attitude maneuver and attitude tweak. Note that the attitude variations, as shown in Figure 8, result in systematic errors 
in the flux time series of target stars, which need to be modeled and removed by the PDC component 10 before the search 
for planets is initiated. 



arcsecond 



Figure 7. The attitude determination algorithm. 



Figure 8. Difference in declination between nominal attitude and attitude solution for second quarter flight data 


4. CONCLUSIONS 


This paper discusses the algorithms of the PPA software component in the science data processing of the Kepler mission. 
The health and performance of the Kepler photometer are assessed by analyzing time series of instrument metrics with 
the track and trend algorithm, and photometer performance reports are generated for each of 84 outputs and the full focal 
plane. With the motion polynomials for each output and each long cadence, the attitude of the Kepler spacecraft is 
determined with high precision at each long cadence from the compressed and smoothed measurements of the centroids 
of images of~16,000 bright target stars across the focal plane. The PPA component is demonstrated to work effectively 
with the Kepler flight data. 

Funding for the Kepler mission is provided by NASA’s Discovery Program Office, Science Mission Directorate. 
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